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Parton showers in Monte Carlo event generators reflect to a certain accuracy our under- 
standing of QCD radiation at all orders. For observables sensitive to interjet energy flow 
in well defined regions of phase space, it has been known for some time that relevant 
all-order dynamics is substantially more complicated than that encoded via angular 
ordering in parton shower algorithms, even to leading logarithmic accuracy. Here we 
investigate the extent of the numerical mismatch between leading logarithmic analyti- 
cal estimates (resummation) and parton showers in an effort to better understand the 
accuracy of parton showers for such observables. 

1 Introduction 

Monte Carlo (MC) event generators are the most important tools at our disposal in the 
search for new physics beyond the standard model. At the large hadron collider (LHC) 
for example, alongside the hard process where the production of hitherto unseen particles 
is expected to occur, there will in general be copious amounts of initial and final state 
QCD radiation, conversion of partons into hadrons (hadronisation), effects due to spectator 
interactions, multiple hard scatterings and pile-up. 

In such a complex hadronic environment as good a modeling as possible of all of the 
above physics aspects, is of clear benefit in order to exploit the full discovery potential of the 
LHC. MC tools exist and are being further developed to provide a simulation of the physics 
described above. However of the whole range of effects mentioned above as accompanying the 
hard process, the only part that we understand directly from first principles is perturbative 
QCD radiation. This understanding is in principle what ought to be reflected in the parton 
shower algorithms in various event generators. 

Another source of information on QCD at all orders is provided by analytical tools such 
as those employed in resummation of large logarithms in QCD. This field too has seen some 
interesting relatively recent developments with the observation that for a large class of ob- 
servables (those measuring flow of hadronic energy into limited regions of phase space and 
hence called non global observables) relevant all-order dynamics is much more complicated 
than previously imagined [2j |3] . This more complex soft gluon dynamics is not fully de- 
scribed by parton showers in event generators such as HERWIG [4] and PYTHIA [5] @. The 
HERWIG parton shower for example is based on angular ordering of soft gluon radiation 
which correctly captures the leading logarithms for several observables. However for non 
global observables the angular ordering approximation is not formally correct at single log- 
arithmic accuracy, which in some important cases is also leading logarithmic accuracy. In 
these cases the HERWIG shower cannot be a priori expected to describe even the leading 
logarithmic perturbative behaviour. This is even more the case for the old PYTHIA shower 
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a This dynamics is in fact included only in the ARIADNE event generator [B] 
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(before v 6.3) where angular ordering is imposed on showers pre-ordered in invariant mass, 
which leads to a poor decription of soft emissions at large angles and consequently non- 
global observables in particular. This is potentially worrying since non-global observables 
such as interjet energy flow are commonly measured experimentally and moreover they are 
used to improve MC models by tuning the model parameters to data on such observables 
e.g. hadronic activity in toward or away regions wrt the leading jet, in hadron collisions. 

Given the prevalence and importance of non-global observables in QCD studies, it is 
worth investigating quantitatively how good a description is provided, in such cases, by 
the different parton showers e.g. those in HERWIG and PYTHIA. This can be done by 
comparing the parton showers to the state-of the-art analytical resummations in an effort to 
see what perturbative accuracy can be relied upon in the corresponding MC descriptions. 

In order to get a clearer idea about what one may expect from parton showers such as 
those ordered in angle we first study to what extent the leading single-logarithmic behaviour 
can be obtained by assuming angular ordering for energy flow into a rapidity interval be- 
tween hard jets. Thus we compare the full single-logarithmic resummation (possible only in 
the large N c limit) with a numerical result obtained by addressing just angular ordered con- 
figurations. We then check if our conclusions are borne out by comparing the resummation 
to the real HERWIG parton shower. We then compare the results to the PYTHIA event 
generator both with the old shower and the new shower and draw some conclusions. 

2 Angular ordering and non global observables 

The observable we shall study is the differential E t flow, ^-j§^ into a rapidity interval of 
width A?7 between a pair of jets. For simplicity we shall take these to be back-to-back 
jets produced in e + e~ annihilation but our conclusions generalise straightforwardly to gaps 
between jets in any hard process. This observable involves other than the scale Et an 
additional scale Q, the hard scale for the process. Typically Et <C Q and one then has to 
deal with logarithms a" ln m Q/E t , with m < n. We note that the leading logarithms in this 
case are single logarithms i.e those with m = n. 

The resummation of the leading logarithms has been possible only in the large JV C limit 
[21 G] . In this limit the complex colour algebra involved in the soft large-angle multi-gluon 
emission probabilities is greatly simplified and can be encoded in a dipole evolution Monte- 
Carlo algorithm, since the emitting ensemble reduces to a system of dipoles. Considering 
a dipole configuration C we can write the probability of emitting a soft gluon at scale 
L' = ]xlQ/lo', where u>' is the gluon energy, as : 



where C is the new configuration of dipoles after the emission of the gluon, A(L', L) is the 
no-emission probability between scales L and L' and F c (9',<fi') is the angular part of the 
dipole real emission probability. Thus explicitly we have 



Then the result that contains the resummation of the leading single-logarithms is a form 



P'ciL') 



a s (L')A c (L',L)Fc(9'(b')Pc(L), 



(1) 
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factor, S, which can be written 



(3) 



In other words the resummed result for the probability that the transverse energy in the 
gap is less than some value E t is obtained by summing over all dipole configurations such 
that the gap region f2 is free from emissions above that scale. In order to study the impact 
of angular ordering we replace F by an angular ordered (AO) approximation q 



F{e k ,4> k ) 



(cos 9^ — cos 8ij ) (cos 9jk — cos 6^ ) 

(l-COS0jfe) (1 -COS^fc) 



(4) 



such that the emission k is produced in cones around the hard legs i and j with angular size 
set by the opening angle of dipole ij. This angular ordering is expected to correspond to 
the piece included in e.g. the HERWIG shower. 

In Fig. [T]we plot the resummed inte- 
grated E t cross-section (equivalently the 
form factor E) as a function of t ~ §Mn -S- 
in both the AO approximation as well as the 
full resummation in the large N c limit, for 
A rj = 1.0. We note that angular ordering a 
has only a modest effect and the results at 
say t — 0.15 corresponding to Q = 100 GeV 
E t = 1 GeV differ by just 10 %. Also shown 
there (the 'primary' curve) is the resumma- 
tion in the independent emission approxi- 
mation assuming soft gluons to be emitted 
by just the hard partons with no correlated 
gluon emission other than effects included in 
the running coupling. The angular ordered 

and full results differ substantially from this stronger assumption. 
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Figure 1: 



3 Comparison to HERWIG and PYTHIA 

We now proceed to compare to actual par- 
ton showers in HERWIG and PYTHIA. „., c 



Since we are interested in the dependence on 
the logarithmic variable t we can in fact go 
to very large Q values, which helps us to kill 
potentially spurious subleading terms of rel- 
ative order a s (Q), that would be present in 
the MC showers. Once we draw conclusions 
for a given value of t , we can apply those 
observations to an appropriate E t value for 




b This replacement would be exactly correct for global obsei'vables where 1 "one ^gijjd fr'flfcly integratc°°8vcr 
azimuthal angles. 

Figure 2: 
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itany Q. Given that angular ordering pro- 
duces results that are in reasonable agree- 
ment with resummation we can a priori ex- 
pect the HERWIG shower to also not differ 
substantially. This is in fact the case as can 

be seen from Fig. [2] where once again for t = 0.15 we note just a 10 % difference between 
HERWIG and resummation (essentially indistinguishable curves on the plot here). 

For the old PYTHIA shower where an- 
gular ordering is imposed on top of order- 
ing in invariant mass the results are sig- 
nificantly different. Here one notes from 
Fig. [3] a discrepancy of 50 % for t = 0.15 
which may be misattributed to other (non- f 
perturbative) physics during the procedure | 
of tuning PYTHIA to data. Bearing in o.c 
mind that the non global logarithmic ef- 
fects arc not universal and are absent for 
other (global) observables this is a potcn- ie 
tially dangerous and inaccurate way to ac- 
count for them. 

For the new PYTHIA shower which im- 
plements angular ordering of soft radiation 

via a dipole type phase-space the results show a marked improvement and a difference of 
only 7.5 % is seen at t — 0.15. However as we increase the width of the rapidity slice in 
question the description from the new PYTHIA shower deteriorates and we observed that 
for At] = 3.0 the new shower performs poorly when compared to resummation and HERWIG 
results. 




Figure 3: 



4 Conclusions 

We can conclude from our studies that while angular ordering as implemented in QCD 
parton showers does not formally give the correct single logarithmic behaviour for non-global 
observables, numerically one does not have to worry too much about the difference between 
angular ordered and resummed results even in this case. However ordering in other variables 
such as invariant mass and angle (PYTHIA old) gives rise to a poor perturbative description. 
This incorrect perturbative estimate may be adjusted for in other non-perturbative aspects 
of the MC model but this is not optimally accurate. We also observed that the new PYTHIA 
shower is substantially better than its earlier counterpart but noted a puzzling discrepancy 
at large values of the interjet rapidity interval Ar] [7], which we hope can be clarified in the 
near future. 

It is imperative in our opinion, especially for observables related to energy flow between 
hard jets, to compare the results from HERWIG and PYTHIA both at the parton and 
hadron level. If large discrepancies are observed at parton level information from analytical 
resummed calculations where available could be exploited to better understand the uncer- 
tainty and limitations of the MC descriptions. 
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